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Abstract. We introduce a model whose thermal conductivity diverges 
in dimension 1 and 2, while it remains finite in dimension 3. We consider 
a system of oscillators perturbed by a stochastic dynamics conserving 
momentum and energy. We compute thermal conductivity via Green- 
Kubo formula. In the harmonic case we compute the current-current 
time correlation function, that decay like i -d//2 in the unpinned case 
and like t~ d / 2 - 1 if a on-site harmonic potential is present. This implies 
a finite conductivity in d > 3 or in pinned cases, and we compute it 
explicitely. For general anharmonic strictly convex interactions we prove 
some upper bounds for the conductivity that behave qualitatively as in 
the harmonic cases. 



1. Introduction 

The mathematical deduction of Fourier's law and heat equation for the 
diffusion of energy from a microscopic Hamiltonian deterministic dynamics 
is one of the major open problem in non-equilibrium statistical mechanics 
[6]. Even the existence of the thermal conductivity defined by Green-Kubo 
formula, is a challenging mathematical problem and it may be infinite in 
some low dimensional cases |13j . Let us consider the problem in a generic 
lattice system where dynamics conserves energy (between other quantities 
like momentum etc.). For x £ denote by £ x (t) the energy of atom x. To 
simplify notations let us consider the 1-dimensional case. Since the dynamics 
conserves the total energy, there exist energy currents j x ,x+i (local functions 
of the coordinates of the system) , such that 

^x(t) = jx-l,x{t) ~ jx,x+l(t) (1) 

Another consequence of the conservation of energy is that there exists a fam- 
ily of stationary equilibrium measures parametrized by temperature value 
T (between other possible parameters). Let us denote by < • >=< • >t 
the expectation of the system starting from this equilibrium measure, and 
assume that parameters are set so that < j x ,y >= (for example if total 
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momentum is fixed to be null). Typically these measures are Gibbs measure 
with sufficiently fast decay of space correlations so that energy has static 
fluctuation that are Gaussian distributed if properly rescaled in space. Let 
us define the space-time correlations of the energy as 

S(x,t) =< £ x (t)S o (0) >- <£ > 2 ■ 

If thermal conductivity is finite, S(x, t) should be solution of the diffusion 
equation (in a proper large space-time scale) and thermal conductivity (TC) 
can be defined as 

By using the energy conservation law JT]), time and space invariance (see 
section [3]), one can rewrite 

ft \ / rt 



<T) = &2^E((/ WOO*) [I Jo,i(s')ds^ 

1 f°° 
= T2E/ .MO)) dt 



(3) 



which is the celebrated Green-Kubo formula for the thermal conductivity 

(cf. HZ]). 

One can see from ([3]) why the problem is so difficult for deterministic 
dynamics: one needs some control of time decay of the current-current cor- 
relations, a difficult problem even for finite dimensional dynamical systems. 
Furthermore in some one-dimensional systems, like Fermi-Pasta-Ulam chain 
of unpinned oscillators, if total momentum is conserved by the dynamics, 
thermal conductivity is expected to be infinite (cf. [13] for a review of nu- 
merical results on this topic) . Very few mathematically rigorous results exist 
for deterministic systems ([15] IB]). 

In this paper we consider stochastic perturbations of a deterministic 
Hamiltonian dynamics on a multidimensional lattice and we study the cor- 
responding thermal conductivity as defined by ([3]). The stochastic pertur- 
bations are such that they exchange momentum between particles with a 
local random mechanism that conserves total energy and total momentum. 

Thermal conductivity of Hamiltonian systems with stochastic dynamical 
perturbations have been studied for harmonic chains. In [5j[7] the stochastic 
perturbation does not conserve energy, and in [3] only energy is conserved. 
The novelty of our work is that our stochastic perturbations conserve also 
momentum, with dramatic consequences in low dimensional systems. In fact 
we prove that for unpinned systems (where also the Hamiltonian dynamics 
conserve momentum, see next section for a precise definition) with harmonic 
interactions, thermal conductivity is infinite in 1 and 2 dimensions, while 
is finite for d > 3 or for pinned systems. Notice that for stochastic per- 
turbations of harmonic systems that do not conserve momentum, thermal 
conductivity is always finite [3] . 
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This divergence of TC in dimension 1 and 2 is expected generically for 
deterministic Hamiltonian non-linear system when unpinned. So TC in our 
model behaves qualitatively like in a deterministic non-linear system, i.e. 
these stochastic interactions reproduce some of the features of the non- 
linear deterministic hamiltonian interactions. Also notice that because of 
the conservation laws, the noise that we introduce is of multiplicative type, 
i.e. intrinsically non- linear (cf. ([6]) and (J7|)). On the other hand, purely 
deterministic harmonic chains (pinned or unpinned and in any dimension) 
have always infinite conductivity [15] . In fact in these linear systems energy 
fluctuations are transported ballistically by waves that do not interact with 
each other. Consequently, in the harmonic case, our noise is entirely respon- 
sable for the finiteness of the TC in dimension 3 and for pinned systems. 
Also in dimension 1 and 2, the divergence of TC for unpinned harmonic 
systems is due to a super diffusion of the energy fluctuations, not to bal- 
listic transport (see [21 [12] where this behavior is explained with a kinetic 
argument). 

For anharmonic systems, even with the stochastic noise we are not able 
to prove the existence of thermal conductivity (finite or infinite). If the 
dimension d is greater than 3 and the system is pinned, we get a uniform 
bound on the finite size system conductivity. For low dimensional pinned 
systems (d = 1,2), we can show the conductivity is finite if the interaction 
potential is quadratic and the pinning is generic. For the unpinned system 
we have to assume that the interaction between nearest-neighbor particles 
is strictly convex and quadratically bounded at infinity. This because we 
need some informations on the spatial decay of correlations in the stationary 
equilibrium measure, that decay slow in unpinned system [9j. In this case, 
we prove the conductivity is finite in dimension d > 3 and we obtain upper 
bounds in the size N of the system of the form in d = 1 and (logiV) 2 
in d = 2 (see Theorem [3] for precise statements) . 

The paper is organized as follows. Section [2] is devoted to the precise 
description of the dynamics. In section [3l we present our results. The 
proofs of the harmonic case are in section 2] and [5] while the proofs of the 
anharmonic case are stated in section [6l The final section contains technical 
lemmas related to equivalence of ensembles. 

Notations : The canonical basis of M. d is noted (ei,e2, ■ ■ ■ , e^) and the 
coordinates of a vector u £ M. d are noted (u 1 , . . . ,u d ). Its Euclidian norm 
|u| is given by |u| = y (u 1 ) 2 + . . . + (u d ) 2 and the scalar product of u and 
v is u • v. 

If N is a positive integer, Z 1 ^ denotes the d-dimensional discrete torus of 
length N and we identify x = x + kNej for any j = 1, . . . , d and k E 7L. 

If F is a function from Z d (or Z^r) into M then the (discrete) gradient of F 
in the direction ej is defined by (V ej i ? )(x) = F(x+ej)— -F(x) and the Lapla- 
cian of F is given by (AF)(x) = Y? i=l {F(x + ej) + F(x - e,-) - 2F(x)}. 
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2. The dynamics 

In order to avoid difficulties with definitions of the dynamics and its sta- 
tionary Gibbs measures, we start with a finite system and we will define 
thermal conductivity through an infinite volume limit procedure (see sec- 
tion Ej). 

We consider the dynamics of the system of length N with periodic bound- 
ary conditions. The atoms are labeled by x G Zjy. Momentum of atom x is 
p x G M. d and its displacement from its equilibrium position is q x G M. d . The 
Hamiltonian is given by 



n N = 



f 4r L 



l ^ + W(c lx )+ 1 - £ F(q x -q y ) 



|y-x|=i 

We assume that V and W have the following form: 

d d 

^(q x - oy) = E V M - 4). w (<t*) = E W M)- 

3=1 3=1 

and that Vj,Wj are smooth and even. We call V the interaction potential, 
and W the pinning potential. The case where W = will be called unpinned. 
We consider the stochastic dynamics generated by the operator 

L = A + jS. (4) 
The operator A is the usual Hamiltonian vector field 
A =E {P* ' 5 <ix - d^JiN ■ <9 Px } 

X 

while 5 is the generator of the stochastic perturbation and 7 > is a 
positive parameter that regulates its strength. The operator S acts only on 
the momentums {p x } and generates a diffusion on the surface of constant 
kinetic energy and constant momentum. This is defined as follows. If d > 2, 
for every nearest neighbor atoms x and z, consider the d — 1 dimensional 
surface of constant kinetic energy and momentum 

Se, P = |(Px,Pz) G M. 2d : i(|p x | 2 + |Pz| 2 ) = e; p x + p z = p| . 
The following vector fields are tangent to § ejP 



so Yli j=i(X^) 2 generates a diffusion on S eiP (see [H]). In d > 2 we define 

I d 
S = 2(d-l) ^2^2{ X ^+ e k 

= 4(d — I) ^ EWz) 

x,z6Z^ 1,3 
|x-z| = l 



THERMAL CONDUCTIVITY 5 

where ei, . . . , is canonical basis of Z d . 

Observe that this noise conserves the total momentum ^ x p x and energy 
Hn, i-e. 

s^ Px = o, sn N = o 

X 

In dimension 1, in order to conserve total momentum and total kinetic 
energy, we have to consider a random exchange of momentum between three 
consecutive atoms (because if d = 1, S ejP has dimension 0), and we define 



6 

where 



Y x = (p x - px+i)dp x ^ + (p x+ i - Px-\)d Px + (pa,-! - Px)d Px+1 

which is vector field tangent to the surface of constant energy and momen- 
tum of the three particles involved. 

The corresponding Fokker-Planck equation for the time evolution of the 
probability distribution P(q, p, t), given an initial distribution -P(q, p, 0) is 
given by 

^ = (-A + 1 S)P = L*P. (5) 

where L* is the the adjoint of L with respect to the Lebesgue measure. 

Let {wl^y] x, y E T^; i,j = 1, . . . ,d; |y — x| = 1} be independent stan- 
dard Wiener processes, such that wl£ y = ro^ x . Equation ([5]) corresponds 
to the law at time t of the solution of the following stochastic differential 
equations: 

dq x = Px dt 

dp x = —O^Hn dt + 27Ap x dt 

d 



+ ^f2=r E ^ E (^ z Px) dw^t) 



(6) 



z:|z-x|=l i,j=l 



In d = 1 these are: 



dp x = -d qx H N dt + -A(Ap x +p x -i +p x +i)dt 
o 

+ Y3 E {Yx+kPx) dw x+ k{t) 
» fc=— 1,0,1 



(7) 



where here {uu x (t),x = 1,...,N} are independent standard Wiener pro- 
cesses. 

Defining the energy of the atom x as 

£x = \vl + ^(q x ) + \ E V ^y - 

y:|y-x|=l 
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the energy conservation law can be read locally as 

d 

£x(i) -£x(0) = J x _ efc ,x([0,t]) " J x ,x + e fc ([M)) 

fc=l 

where J x ,x+e fc ([0, £]) is the total energy current between x and x + up to 
time t. This can be written as 

Jx,x+e fc ([0,*]) = / ix,x+e fc (s)d S + M XiX+efc (i) (8) 
J 

In the above M XiX+efe (t) are martingales that can be written explicitly as 
Ito stochastic integrals 

M X)X+e , {t) = £ jf * (^ x f x+e /x) (*) ^< x+ e fc (-) (9) 

In ci = 1 these martingales write explicitly as 

M X)X+ i(t) = * H- Y] (Y x+k £ x ) dw x+k (t) 

V 3 J ° fc=-l,0,l 

The instantaneous energy currents j XiX + efc satisfy the equation 



(10) 



Jx,x+e fe = -o( VF )( f lx+e fc - qx) • (Px+e fc + P> 



d 

L£ x = (jx-e k ,x — Jx,x+e fc ) 
fc=l 

and it can be written as 

jx,x+e k = 3x,x+e k + 7Jx,x+e fe ■ (H) 

The first term in (jllj) is the Hamiltonian contribution to the energy current 

^< 

d (12) 
j'=i 

while the noise contribution in d > 2 is 

7ix,x+e fc = -7(V efeP 2 )x (13) 

and in d = 1 is 

lfx,x+l = - "fV(f(Px-l,Px,Px+l) 
<p(Px-l,Px,Px+l) = \\PI+1 + ^Pl +PI-1 +Px+lPx-l ~ 2Px+lPx ~ 2p x p x -l\ 

In the unpinned case (W = 0), given any values of £ > 0, the uniform 
probability measure on the constant energy-momentum shell 



SjV,£ = < (P, q) : n N = N£, = 0, q x = 



xei 
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is stationary for the dynamics, and A and S are respectively antisymmetric 
and symmetric with respect to this measure. For the stochastic dynamics, 
we believe that these measures are also ergodic, i.e. total energy, total mo- 
mentum and center of mass are the only conserved quantities. Notice that, 
because of the periodic boundary conditions, no other conserved quantities 
associated to the distortion of the lattice exist. For example in d = 1 the 
total length of the chain ^2 x (q x +i — Qx) is automatically null. 

In the pinned case, total momentum is not conserved, and the ergodic 
stationary measures are given by the uniform probability measures on the 
energy shells 

S JV)£ = {(p,q) :H N = NS}. 

In both cases we refer to these measures as microcanonical Gibbs mea- 
sures. We denote by < • > N ,e the expectation with respect to these micro- 
canonical measures. 

We will also consider the dynamics starting from the canonical Gibbs 
measure < • >jv iT with temperature T > defined on the phase space 
(M. 2d f d N by 

e -Wjv/T 
< • > n ,t= — y dqdp. 

To avoid confusions between these measures we restrict the use of subscript 
£ for the microcanonical measure and subscript T for the canonical measure. 

3. Green-Kubo formula and statement of the results 

In the physical literature several variations of the Green-Kubo formula ([3]) 
can be found ([l3j[7])- As in ([3]), one can start with the infinite system and 
sum over all x 6 Z . One can also start working with the finite system with 
periodic boundary conditions and sum over x £ where A^ is a finite 
box of size iV and take the thermodynamic limit N —* oo (before sending 
the time to infinity). In the finite case there is a choice of the equilibrium 
measure. If < • > is the canonical measure at temperature T, one refers 
to the derivation a la Kubo. If < • > is the microcanonical measure at 
energy £N d , one refers to the derivation a la Green. Because of to the 
equivalence of ensembles one expects that these different definitions give 
all the same value of the conductivity, provided that temperature T and 
energy £ are suitably related by the corresponding thermodynamis relation. 
Nevertheless a rigorous justification is absent in the literature. 

In the sequel we will consider the microcanonical Green-Kubo formula 
(noted k) and the canonical Green-Kubo formula (noted k) starting from 
our finite system. 

In the harmonic case we work out the microcanonical Green-Kubo version 
that we compute explicitly. Similar computations are valid (with less work) 
for the canonical version of the Green-Kubo formula and will give the same 
result. In the anharmonic case equivalence of ensembles is less developed 
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and we deal only with the canonical version of the Green-Kubo formula. 

The microcanonical Green-Kubo formula for the conductivity in the di- 
rection ei is defined as the limit (when it exists) 

k^CT) = lim lim -L £ E„, £ [Jx, x+ei ([0, t]) J Q , ei ([0, t])] (14) 

where E Nj£ is the expectation starting with the microcanonical distribution 
< • >jv i£ , and the energy £ = £{T) is chosen such that it corresponds to 
the thermodynamic energy at temperature T (i.e. the average of the kinetic 
energy in the canonical measure). In the harmonic case T = £. 

Similarly the canonical version of the Green-Kubo formula is given by 

R 1 > 1 (T)= lim lim -L E^ T [J x , x+ei ([0,t])J , ei ([0,t])] (15) 



when this limit exists. Here M NT indicates the expectation with respect to 
the equilibrium dynamics starting with the canonical measure < • >jv jT at 
temperature T. These definitions are consistent with ©-([2D as we show at 
the end of this section. 

Our first results concern the (a, ^-harmonic case: 

Vj(r)=ar 2 , Wj(q) = uq 2 , a > 0, v > (16) 

Theorem 1. In the (a, v) -harmonic case H6]) . the limits defining k ' and 
K ' exist. They are finite if d > 3 or if the on-site harmonic potential is 
present (v > 0), and are infinite in the other cases. When finite, k(T) and 
k(T) are independent ofT, coincide and the following formula holds 



k l ^(T) = k ' (T) = — =— / v , ' v ' dk + -L (17) 
8ir*d>yJ m * V(k) d 1 ; 



where cj(k) is the dispertion relation 



w (k) = [ v + 4a sin 2 (7rk- 7 ') J (18) 



and 



(k)= ,8El =1 sin 2 (^), i/ d>2 

M/3 sin 2 (vrk)(l + 2cos 2 (7rk)), i/ d = l 

Consequently in the unpinned harmonic cases in dimension d = 1 and 2, 
the conductivity of our model diverges. In order to understand the nature 
of this divergence we define the (microcanonical) conductivity of the finite 
system of size N as 



1 1 



■/x.x+exdMjv]) 



7d 



(20) 
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where ijv = N/v s with v s = lim^o (c^o^k)! = 2a 1 / 2 the sound velocity. 
This definition of the conductivity of the finite system is motivated by the 
following consideration: Vfctj(k) is the group velocity of the k-mode waves, 
and typically v s is an upper bound for these velocities. Consequently t]y is 
the typical time a low k (acoustic) mode takes to cross around the system 
once. One defines similarly (T) by 



1 1 

xe 



^ <^x,x+ei([0,ijv]) 



(21) 



We conjecture that kn (resp. Rn) has the same asymptotic behavior 
as the conductivity defined in the non-equilibrium stationary state on the 
open system with thermostats at the boundary at different temperature, as 
defined in eg. El [15] . 

With these definitions we have the following theorem: 

Theorem 2. In the harmonic case, ifW = 0: 

(1) k n ~ N 1 / 2 ifd = l, 

(2) k n ~ logiV ifd = 2. 

In all other cases kn is bounded in N and converges to k. 
Same results are valid for kn. 

In fact we show that, in the harmonic case, we have 

lim^p^Ll (22) 

This is a consequence of equation (|5ip that one can easily check is also valid 
if the microcanonical measure is replaced by the canonical measure. 

In the anharmonic case we cannot prove the existence of neither k ' (T) 
nor k ' (£ ), but we can establish upper bounds for the canonical version of 
the finite size Green-Kubo formula (|2ip . Extra assumptions on the poten- 
tials V and W assuring a uniform control on the canonical static correlations 
(see ()86H89p ) have to be done. In the unpinned case W = 0, ([89]) is valid as 
soon as V is strictly convex. In the pinned case W > 0, (|86p is "morally" 
valid as soon as the infinite volume Gibbs measure is unique. Exact assump- 
tions are given in [3], theorem 3.1 and theorem 3.2. In the sequel, "general 
anharmonic case" will refer to potentials V and W such that ([861) (or ([891) ) 
is valid. 

Theorem 3. Consider the general anharmonic case. There exists a constant 
C (depending on the temperature T) such that 

• Ford> 3, 

(1) either W > is general 

(2) orifW = and < c_ < V" < C + < oo for any j, 
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then 

~4\T)<C. 

• For d = 2, if W = and < c_ < V" < C+ < oo for any j, then 

~4\T)<C(logN) 2 . 

• For d = 1, if W = and < c_ < V" < C+ < oo, then 

K^(T)<CVN. 

• Moreover, in any dimension, if Vj are quadratic and W > is gen- 
eral then KjX (T) < C. 

The proof of this statement is in section [6J 

We now relate the definition of the Green-Kubo (j!4j) and (|15p to the vari- 
ance of the energy-energy correlations function ^j. 

Consider the infinite volume dynamics on % d under the infinite volume 
canonical Gibbs measure with temperature T > 0. The expectation is de- 
noted by E T . Fix t > and assume that the following sum makes sense 

D${t) = Y, * 1 ^t [ - T)(£b(0) -T)]=J2 x i x^(x, t). (23) 

xGZ d xSZ d 

If x 7^ 0, by space and time invariance of the dynamics, we have 

E T [(£ x (i) - T)(f (0) - T)] = -^E T [ (£ x (i) - £ x (0)) (£b(t) - £ (0)) ] (24) 

By definition of the current, we have for any y E Z d : 

d 

£ y (t)-£ y (0) =^(J y _ efe , y ([0,i]) - J y , y+efe ([0,i])) (25) 

k=l 

By two discrete integration by parts one obtains 

D?(t)=Y E T [ Jx,x+e, ([0, t]) J ,e 3 ([0, *] )] (26) 

a;eZ d 

so that the thermal conductivity is equal to the space-time correlations of 
the total current 

k^(T) = 5 (i - j) km ^ E E - [ J x,x +ei ([0,t]) Jo,e,([0,t])] (27) 

x&L d 

Of course this derivation is only formal even for fixed time t > 0. The 
problem is to define the infinite volume dynamics and to show S(x,t) has 
a sufficiently fast decay in x. For the purely Hamiltonian dynamics, it is 
a challenging problem. For the stochastic dynamics it seems less difficult 
but remains technical. To avoid tkese difficulties we adopt a finite volume 
limit procedure starting from ([3]). Tkis explains tke definitions (jl4h and (|15p . 
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Consider now the closed dynamics on 7L N starting from the microcanonical 
state. The rest of the section is devoted to the proof of the following formula 



and an identical formula in the canonical case (with E Nj£ substituted by 



The term 7/din (|28h is the direct contribution of the stochastic dynamics 
to the thermal conductivity. In the microcaconical case we actually prove 
that is is equal to j/d only for the harmonic case. A complete proof of 
(|28p for anharmonic interaction demands an extension of the equivalence 
of ensembles estimates proven in section [7J In the grancanonical case this 
problem does not appear. 

Starting in the microcanonical case, remark that the first term on the 
RHS of (|28p can be written as 



If 7 = 0, which corresponds to the purely Hamiltonian system, as N 
and then t goes to infinity, and if one can prove that the current-current 
correlation function has a sufficiently fast decay, then one recovers the usual 
Green-Kubo formula (J3]). 




(28) 



E NjT ). 




(29) 
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To prove (|29l) one uses space and time translation invariance of the dy- 
namics 

2 N 

/ r I 

-i2j. 7\rd\ — It 



(2T z tN d )- 1 E N ,£ 



Jx.x+ex ( s )ds 

= (2ThN d Y 1 Y, f ds f du Ejv,£ (j^. x+ei (*) iy, y+ei («)) 

x,y -JO JO 

= {T 2 tN d r l Y I ds f du E Ni£ (i« x+ei (s) % >y+ei («)) 

= {T 2 tN d )~ l Y I ds f E^, £ (i5_ y , x _ y+ei (a - u) ioV (0)) 

= ^jf (i--J Y E ^(^ +ei isK, e M ds 

We now give the proof of (|28|) . Because of the periodic boundary condi- 
tions, since j s if a gradient (cf. (|13p ). the corresponding terms cancel, and 
we can write 



Ex^x,x +ei ([0,t]) = /o t Exi5 I x+e 1 (*)^ + Ex^x,x+e 1 (t) 

= Jo3 ei W ds + ^ ei (t) 

so that 

2 N 



(30) 



^Jx,x +ei ([0,t]) 




n 2^ 



3ei(s)rfs 



+ (jtN d )- 1 M Nte {m 2 ei (t)) 



(31) 



3d (s) ds 



The third term on the RHS of (|3ip is shown to be zero by a time reversal 
argument and the second term on the RHS of pip gives in the limit a 
contribution equal to j/d. 

To see the first claim let us denote by {w(s)}o< s <t the process 
{(p x (s), q x (s)); x G Zfj, < s < t} arising in (|6|) or in ([7]) for the one- 
dimensional case. The reversed process {uj*}o< s < t is defined as to* = 0Jt- s - 
Under the microcanonical measure, the time reversed process is still Markov 
with generator —^4 + 75. The total current Jt(w.) = E x ^x,x+ei ([0, t]) is 
a functional of {w s }o< s <i. By ([6HTJ) , we have in fact that Jt(-) is an anti- 
symmetric functional of {to s }o< s <t, meaning 



Jt(Ws}o<s<t) — —Jt(Ws}o<s<t) 



(32) 
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In fact, similarly to (J5]), we have 

J s (lo*)= I'* ti ei )*(LU*(v))dv + Wl* ei (s), 0<s<t (33) 

J 

where (SD r l* 1 (s))o<s<t is a martingale with respect to the natural filtration 
of (uJ* s )o< s <t and (3 e J* = Ex(i a )i,x+ ei is equal to -3 ei = -Exi5,x+ei- 
We have then by time reversal 

E NtE [J t (u.)3 ei (u;(t))} = -E N , £ [J t (u;*)3 ei (u;*(0))} 

= -E N , f U (3 ei )* (u* («))ds + (t)^ 3 ei (u* (0)) 

(34) 

where the last equality follows from the martingale property of 9JT . Recall 
now that (3ei)* = — 3ei- By variables change s — ► f — s in the time integral, 
we get 

ft 



3^*(w*(s))efo 3eiK(0)) 



E JVi£ [J t (o;.)5e 1 (a;(t))] = E JV , £ 
It follows that 

3 ei (w(s))(is^ an ei (i) 



3 ei (w(s))ds 3ei(w(t)) 



E 



JV.£ 







E 



5ei(w(s))9Jtei(s)ds 







iV.f 







3 ei (w(s)) ( J s (w.) - y 3 ei (w(v))dv 
For the second term on the RHS of (1311) we have 

r<^— lie /'<nrf2 _ 7 

lK}) (d-l)N d - 

x 1,3 



(35) 



(36) 



( i iV d )- 1 E,, £ «(t)) = (d _[ )Nd EE{ (^fx +ei (Px/2) 



(d 


- l)JV d 




7 


(d 


- l)JV d 




2 7 


(d 


- l)jV d 




27 



n N d ^ ^ ( (^*+ e i PxPx+ei 



l)]yd ^ ((PxPx+e, PxPx+ ei , 



N,£ 
2 



(d~ l)N d — \ /iV,£ 

Thanks to the equivalence of ensembles (cf. lemma [7]), this last quantity is 
equal to 

27 T + N ' d ° N (37) 

where On remains bounded as N — ► oo. The calculation in d = 1 is similar. 
The contribution of the martingale term for the conductivity is hence 7/d 
and we have shown (f28l) . Notice this is the only point where we have used 
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the equivalence of ensembles results of section [7] that we have proven only 
in the harmonic case. We conjecture these are true also for the anharmonic 
cases. 



Observe that all the arguments above between (|30j) and (|37j) apply di- 
rectly also to the canonical definition of the Green-Kubo but without the 
small error in N (because for the canonical measure momentums p x are 
independently distributed and the equivalence of ensembles approximations 
are in fact equalities). Therefore we have the similar formula to (128j) : 



^ Yl E ^ (</x,x +ei ( [0, t] ) J„,ei ( [0, t] ) ) 

Y I &,*+*M ds 

Z Jo 



2TH 

X 

= (2T 2 N d t)' l E Ni 



+ 



1 



(38) 



In the next sections we will consider the (a, z^)-harmonic case and we will 
compute explicitly the limit (as N — > oo and then t — ► oo) of the two first 
term on the RHS of (15TD. 



4. Correlation function of the energy current in the 

harmonic case 

We consider the (a, ^-harmonic case (fl~6j) . We recall that 3ei = X)x.?x,x+ei • 
Because of the periodic boundary conditions, and being j£ x +ei a spatial gra- 
dient (cf. (PT3j) ). we have that $ ei = Ex^xx+ei- We are interested in the 
decay of the correlation function: 

C ltl (t) = lim -L Ejv , £ (3 ei (^ ei (0)) = lim ^E^ £ (j %(0)j° x+ei (i)) 

(39) 

where Ejy |£ is the expectation starting with the microcanonical distribution 
defined above. 

For A > 0, let ma,tv be the solution of the Poisson equation 



\U X ,N ~ Lu X ,N = -X^x,x+ei 



given explicitely in lemma [2] of section [5j By lemma [H we can write the 
Laplace transform of Cn(f) as 



tfte-^Ci,!^) (it = lim (j$, ei ux, N ) N£ (40) 
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Substituting in (flOl) the explicit form of ua,at given in lemma El we have: 
Op 1 _ — ^ 

- \3o, ei u X,N} N>£ = ^-2^5A,Af(x-y) ((q ei - q ) • (p 6l +Po)(Px ■ Qy)) N>£ 

a 2 

= ^"X)ffA,iv(x - y) ((q ei • po - q • PeJ(Px ■ <ly)) N , e 

Oi 2 

+ 2^X)fA,iv(x-y) ((q ei -p ei - q • p )(Px • 

x,y 

(41) 

Observe that the last term on the RHS of (|4ip is null by the translation 
invariance property. So we have (using again the translation invariance and 
the antisymmetry of <?a,jv) 

a 2 

~ (io,e 1 «A,Jv) JVe = 2Z%29\n{x ~ y) ((qei - q-ei) ■ Po)(Px ■ 
x,y 



Define 



2- 

X 



In the unpinned case v = 0, conditionally to the positions configuration q, 

the law of p is // q = fi A 

X 

measure on the surface 



the law of p is /x q = fi N , — (defined in lemma [6]), meaning the uniform 



IPxJxez^; 



\^2pI = *0v(q); ^ Px = f 

X X J 



By using properties (i),(ii) and (iii) of lemma [6l one has for 

(((qei - qej • po)(px • qy))^ = Yl (poPx) (qli - q-ei)q y > Ni£ 

hi 

= E (popi) (< - i-=i)4>„, e 

' / ' 2K NM V («) 

^\dA^(iV rf -l) iqei 4 " ei;4y / iV , £ 

1 d 

= -]yd3iE<(Po) 2 (qei -q-ejqy)^ 

8=1 
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For x = 0, one gets 

<((q ei - q ei ) • po)(po • <hr)) N>e = Yl {f** (poPo) (qli - q-ejqy 



i=l 



<^ q (p^,po) « - q- ei )q^) JVi£ (43) 



i=l 



iV.fi 



In the pinned case ^ > 0, conditionally to the positions configuration q, 
the law of p is A q = X N , — — — (defined in lemma [5]), meaning the uniform 

V 2A'jv (q) 

measure on the surface 



(Px), 



^Epx = ^(q)| 



We proceed in a similar way and we observe that if x ^ 0, Aq(pQp^) = 0.(cf. 
ii) of lemma [5]) 

Since is antisymmetric (see (|64H65p ) and such that Y1z9*,n( z ) = 0, 
one obtains easily in both cases (pinned and unpinned) 



|-E^(y)E((Po) 2 (4-^ ei K 

y i 



N, 



iv.fi 



+ 



2 7 j\T d _ 



I 9\Ay- x )J2(( p °) 2 ( q ei -9- ei )4 



li/=o \ a 2 



JV,£ 



\ 2 



Let rjv(x), x G Z^-, be the unique solution of 

(yl - aA) Tat = <5 ei - £_ ei 

such that X^xez^ F./v(x) = 0. 

By (iii) of lemma [7J and ([77]) , we have 



, l„=o \ a 2 £ 2 



^2g\,N(y)^N(y) 



„ ClogiV ^, „ ClogiV /^V 2 ^- 

-■■ 2^lgA,Jv(y)| < ^ ^(SA.aK*)) j < 



(44) 



(45) 



N d 



a log iv 



(46) 



Hence the last term of (|4"6j) goes to 0. 
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Taking the limit as N — > oo we obtain (see ([80]) ) 

I e- M C 1 , 1 (t)dt = —'£g x (z)T(z) (47) 

where g\ are solutions of the same equations as g\ : n but on Z d and T is the 
solution of the same equation as Tjy but on Z d . 

Using Parseval relation and the explicit form of the Fourier transform of 
9\ (cf. (I74j) ) and T, one gets the following formula for the Laplace transform 
of Ci,i(t) for d > 2: 

a 2 £ 2 /■ / sin 2 (27rk x ) \ 1 

— r- / rfk J 5 s s 48 

d ./[o,l]* V^ + 4aE J= ism 2 (7rki)/ A + 8 7 E J= isin 2 (7rki) 

By injectivity of Laplace tranform, Ci i(i) is given by: 



, a 2 £ 2 /■ / sin 2 (27rk 1 ) \ 

Cn(t) = — — / dk -J J exp< 

d i[o,i] d \^ + 4aX;j = isin 2 (7rkJ)y 



-87t^sin 2 (vrk J ') 
(49) 

For the one dimensional case, the equation for (resp. g\) is different 
(see (I75p ) and we get the following integral representation of the correlation 
function of the energy current: 

C 1A (t) = a£2 dA;cos 2 (7r/c)exp|-^sin 2 (7rA;)(l + 2cos 2 (7rA;)| (50) 



Cx,i(t) = ^nf (d kl u,(k))V^ k )dk (si) 



In any dimension, we have the following unified formula for C\ \{b) 

£ 2 _ 

' '[0,1] 

where w(k) is defined by (fT8|) and ^>(k) by (p~9|) . Observe that the same 
formula holds if we replace K Ni£ by E W?T . In this last case, the situation is 
simpler since we do not need equivalence of ensembles. 

Standard analysis shows the behavior of Ci,i(i) as t goes to infinity is 
governed by the behavior of the function (<9 k iu;(k)) 2 and ip(k) around the 
minimal value of tp which is 0. In fact, ?/>(k) = if and only if k = or 
k = (1, . . . , 1). By symmetry, we can treat only the case k = 0. Around 
k = 0, ?/>(k) ~ a|k| 2 and (<9 re iw(k)) 2 ~ b(y + |k| 2 ) _1 (k 1 ) 2 where a and b are 
positive constants depending on v and a. Essentially, C\ t i(t) has the same 
behavior as 

, (A .l)2 e -«7t[k[» ! , (k l)2 e -a 7 |k|2 

/ dk^ = — — / dk ^ ' —ft 52 

Ae[o,i] d f + |k| 2 td/2+1 J[o,Vt^ ^ + *- 1 |k| 2 
Hence, we have proved the following theorem 
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Theorem 4. In the (a, v) -harmonic case, the current- current time correla- 
tion function Ci } i(t) decays like 

• Ci,i(i) ~ t~ d l 2 in the unpinned case (v = 0) 

• Cn(t) ~ i~ d / 2 - 1 j n ^ e pi nne d case (u > 0) 

5. Conductivity in the harmonic case 

Lemma 1. Consider the (a, v)-harmonic case. For any time t, the following 
limit exists 

Ci,i(t) = lim * E^ £ (3 ei (t)a ei (0)) (53) 

N^oo iV" 

and 

dte- xt C hl (t) dt = lim (jS, ei u x , N ) N£ (54) 
same result holds with Ejv,£ replaced by Ejv,t- 
Proof. We only prove this lemma in the microcanonical setting. Let us define 

/iv(i) = ^E J v, £ (5e 1 (^e 1 (0)) (55) 

We first prove the sequence (/iv)^ is uniformly bounded. By Cauchy- 
Schwarz and stationarity, we have 



I 

Jo 



IM*)I < W^flM N ,e^(^)) N ,s (56) 

We now use symmetry properties of the microcanonical ensemble to show 
this last term is bounded above by a constant independent of N. 

N d < 3e x >N,e = ^2 < io,eiix,x+ei >N,e 
x 

2 d 

= ^ ((<4 - qo)(qx +ei - qi)(Pei + Po)(Px +ei + pL) ns 

x i,j=l 

In the unpinned case v = 0, conditionally to the positions configuration q, 
the law of p is fin = /j, N t — (defined in lemma [6]). 

^/2Xjv(q) 

By using properties (i),(ii) and (iii) of lemma [HJ one has 



2 " 

N- d < Z 2 ei > N ,e= (« " qo)(qk " q- ei - 3q^ + Sq^pj,) 2 }^ £ 

i=l 

(57) 

By Cauchy-Schwarz inequality, the modulus of this last quantity is bounded 
above by 

a[8 < £l > N , £ + l - < £l > N ,e] = y < So >n,e (58) 

where the last equality is a consequence of the invariance by translation of 
< • >at,£. Let (Xi, . . . ,X N d) be a random vector with law A^^, meaning 
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the uniform measure on the iV^-dimensional sphere of radius V N d £. The 
vector of energies (£ x > x G Z, d N ) has the same law as (Xf,...,X"^ d ). By 
lemma HI K(X^) =< £q > jv £ is bounded above by a constant independant 
of N. Hence there exists a positive constant C such that 

\f N (t)\<C (59) 

Similarly, inequality ([59]) can be proved in the pinned case v > 0. Let f(t) 
be any limit point of the sequence (/jv(i))jv>i and choose a subsequence 
(Nk)k>o such that (/jv fc ) converges to / (for the pointwise convergence topol- 
ogy). By Lebesgue's theorem, we have 

POO POO 

lim / e - xt f Nk (t)dt= / e~ xt f{t)dt (60) 
k ^°° Jo Jo 

But we have that 

oo 

e~ Xt f N (t)dt = - < j , Bl , u X:N >N,e (61) 



and we have seen in section [J] this last quantity converges as N goes to 
infinity to 

oo 

e^foo^dt (62) 
where foo is given by (see (I18H19P for the notations) 

/«>(*) = -TT-A (d^MVfe-^dk (63) 

By injectivity of the Laplace transform, we get f(t) = fao(t). Uniqueness of 
limit points implies (/at(*)).n>i converges to foo{t) for any t. It follows also 
we can inverse time integral and infinite volume limit in the left hand side 
of (|54p and the lemma is proved. □ 



Lemma 2. (Resolvent equation) 

U X ,N = (A - L)' 1 I - ix, X +ej I = - XI f A,JV(X - y)Px ■ 
V x / x,y 

where g\ t N{ z ) is the solution (such that ^2 z g\,N( z ) = 0) °f the equation 
2A 

— 5AiiV (z) - 4A 5Aj7V (z) = (5(z + ei) - 5(z - e x )) (64) 

7 

for d > 2, or 
2A 1 

— 9\M Z )- o A [ 4 5A,7V(^) + gA,Jv(* + 1) + ffA,Jv(* - 1)] = (S(z+1)-S(z-1)) 
7 a 

(65) 

/or d = 1. Moreover, Au\^ = and Lu\ } n = lSu\ t N- 
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Proof. We only give the proof for the dimension d > 2 since the proof for 
the one dimensional case is similar. Let u\ : n = ^ Yl x y 9x,n(*- — y)Px ■ qy 
The generator L is equal to the sum of the Liouville operator A and of the 
noise operator 7s*. The action of A on u\ : n is null. Indeed, we have: 

Au x ,n = -^[(aA-!//)qx]' I ^5a,at(x - y)q y J H — ^ 5A,Af(x-y)p x -p y 

^ x \ y / ^ y,x 

(66) 

Here, and in the sequel of the proof, sums indexed by x, y, z are indexed by 
Zjy and sums indexed by i,j, k, £ are indexed by {1, ... , d}. Summation by 
parts can be performed (without outcoming boundary terms since we are 
on the torus) and we get 

Au x , N = -^[(° A - ^.wKx- y)qxq y + - ^3a,jv(x - y)p x - p y (67) 

^ x ^ y,x 

Remark now that the function 5(- — e±) — 5(- + ei) is antisymmetric. Hence 
g\ : N, and consequently Ag\ t N, is still antisymmetric. We have therefore 
Au\^n which is of the form: 

Au x ,n = £{ai(x - y)p x • p y + a 2 (x - y)q x • q y } (68) 
x,y 

with ai,a2 antisymmetric. Using the antisymmetricity of a\ and 02, it is 
easy to show that the last two sums are zero and hence Au\ : n = 0. 

A simple computation shows that if I £ {1, . . . , d} then 

*(px) = 27^Tj£ £(^,WVx) 

- 2 ( X i>e \ 2 (n^ I 2 VfrW ] 2 (r> e ) 

2(d — 1) V^ l x,x+e fc ; VPx; 2(^ — 1) v vi x-e fc ,x; IPx, 1 

= ^Tl E {(Px+e t " Px) " (Px " Px-e fc )} 

= 2A(p4) 

Since the action of S is only on the p's, we have 

lSu x , N = a£5A,Jv(x-y)S , (p x ) -q y 
x,y 

= 2a z ^ 5 r AiA r(x-y)(Ap x ) • q y 
x,y 

= 2a z ^(A 5AiA r)(x-y)p x -q y 
x,y 
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where in the last line, we performed a summation by parts. Since gx,N is 
solution of (f64|) . we have 

Au A ,AT - jSux,N = ^ Px ' (°l x + e i ~ ^-ei ) = ~ E Jx.x+ei (69) 



2 

x 



□ 



Let us define the Fourier transform v(£), £ £ Z^, of the function v on Z"^ 

as 

= v(z) eocp(2<7rf • z/JV) (70) 

The inverse transform is given by 

v ^ = W E "(0 exp(-2i< • z/iV) (71) 

On Z rf we define similarly: 

(k) = v( z ) exp(2ivrk • z), k € [0, l} d (72) 



z 

and its inverse by 



i>(z) = / u(k) exp(— 2i7rk • z) 

7[o,ii d 



For A > 0, the function g x : Z d -> R is the solution on Z d of the equation 
2A 

— g x (z) - 4Ag x (z) = 8 (z + e x - <5 (z + ei), d > 2 (73) 

7 

2A 1 

— 5A (z) - -A (4g x (z) + g x (z + 1) + 5a (z - 1)) = 5 (z + 1) - <5 (z - 1) 
7 3 

(2= 1 

Then we have 

« « x — 2i7rsin(2-7rk 1 ) „ , 

5A(k) = ^ > , if rf> 2 (74) 

— +16EliSin 2 (vrkJ) 

and 

0a(*O = — — , ifd=l (75) 

h -sin 2 (vrA;) (1 + 2 cos 2 (vrA;)) 

7 3 

Since gx,N is the solution of the same equation as gx but on Z 1 ^, we have 
the following formula for gx,N- 

9x,n(0 = 9x(Z/N) (76) 
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The following bound follows easily from Parseval relation: 

E(^(x)) 2 



* < 7 



A 2 



(77) 



Similarly, the function T^v defined in (|45|) has Fourier transform given by 

f n(0 = t(£/N) (78) 

where 



f(k) 



-2i sin^vrk 1 ; 



(79) 



^ + 4aE J= iSin 2 (7rkJ) 
Let us denote by z* the conjugate of the complex number z and observe 



that the function k € [0, l] d g A (k) f (k) 



G M + is continuous. Hence we 



have the following convergence of Riemann sums 



(80) 



Af^oo 



[0,l] d 



dkfc(k)(f (k)]* = £ 5A(y)r(y) 



y e2 



The limits as A — > of the above expressions give the values for the 
conductivity (up to a multiplicative constant) when this is finite. If v = 
it diverges if d = 1 or 2. 



6. ANHARMONIC CASE: BOUNDS ON THE THERMAL CONDUCTIVITY 

We consider in this section the general anharmonic case and we prove 
theorem [3j Recall (|38p . then all we need to estimate is 



(2T Z N 



2^+i r i Ejvr 



V] / ix,x+ ei ( s )^ s 



(81) 



Let us define X^xixx+ei = 3en then we have the general bound (|16j. 
lemma 3.9) 



E 



N,T 



< 10N(3 ei ,(N- 1 - 7 S')- 1 5e i ; 



(82) 



1 



Recall that 5(p x ) = 2A(p x ) if d > 2 and 5(p x .) = -A(Ap x +p x+ i +p x -i) 

6 

if d = 1. 



(iV^-T^)- 1 ^ =EE G ^( X -y)P^'K+e 1 -q'y) (83) 

j=i y 
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where Gjv(z) is the solution of the resolvent equation 

N- 1 G N (z)-2 1 (AG N )(z) = ~[5 (z)+5 ei (z)], d>2 

N^GNiz) - I [A(AG N )(z) + (AGW)(z + 1) + (AG N )(z - 1)] 
o 

= -i[5 ( z ) + ^(z)], d=i 

v Z 

(84) 

The left hand side of (|82p is equal to 

-5TiV^^;^(G J v(x) + G J v(x + e 1 ))(^Vx + e 1 -^(<-^ Q )) JVT 

(85) 

• Pinned case 

In the pinned case, the correlations ( V'(q J x+e — qi.)V'(q J ei — ql) ) 
decay exponentially in x 

< Ce" c|x| (86) 

It follows that the previous expression is bounded by 
CTW^|G N (x)+G N (x + ei )|e- c l x l 

X 

Since Gn is bounded in d > 3, it follows that (fSTj) is uniformly 
bounded in N. In low dimensions, our estimate are to rough and we 
obtain only diverging upper-bounds. Nevertheless, if Vj(r) = otjr 2 
are quadratics and Wj are general but strictly positive then 

\ (87) 

= aj 1 2 < q£q£ >^ T - < q^_ ei q^ >^, T - < q^ +ei q^ > N , T j 

As a function of x, this quantity is a Laplacian in the first direc- 
tion and by integration by parts, the left-hand side of (|8Tj) is upper 
bounded by 

c£ l(AGW)(x) + (AGW)(x + ei)| e~ c|x| (88) 

X 

By lemma El this quantity is uniformly bounded in N . 

• Unpinned case 

In the unpinned case, we assume that < c < V" (q) < C <+oo. 
We have (cf. [9], theorem 6.2, that can be proved in finite volume 



(V J , (^ + e 1 -9x)V J '«-^)) 



24 GIADA BASILE, CEDRIC BERNARDIN, AND STEFANO OLLA 

uniformly) 



^(4 +ei -4)Wei-<?0), 

AT,T 



< C|x|- d (89) 



In the one dimensional case, the random variables r x = q x +i — q x are 
i.i.d. and < V'{r x ) > N , T = 0. Only the term corresponding to x = 
remains in the sum of (f85j) . By lemma El we get the upper bound 

(GW(0) + Gjv(l)) {V'{rl)) NT < C^/N (90) 

For the unpinned two dimensional case, we obtain the upper 
bound 

-d 



C |GW(x) + GW(x + ei)||x| 



< ClogiV \ x \~ 2 ( 91 ) 



~ C(logiV) 2 

For the case d > 3, we use the first point of lemma El ([89]) and 
the fact that 

Y |x|- d ~logiV (92) 

Lemma 3. Let Gm be the solution of the discrete equation {8$ - There 
exists a constant C > independent of N such that 

• Gjv(x) < Cdxl^ 2 + N- 1 ' 2 ), d>3 

• Gjv(x) < ClogiV, d = 2 

• Gjv(x) < CVN, d = l 

• \G N (x + ei) + Gjv(x - ei) - 2G N (x)\ <C, d>l 

Proof. In the proof, C is a constant independent of N but which can change 
from line to line. We first treat the case d > 3. We use Fourier's transform 
representation of Gn: 

i „-2i7rk-x/iV 



where 9n(u) = N 1 + 87 X^j=i sin 2 (7ru- ? ). Gat can also be written in the 
following form 

G N (x) = - ~ [Fjv (x) + F N (x - ei)] (94) 

where 

^ -2i7rk-x/A r 

*W = MEW ,%) 
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Let us introduce the continuous Fourier's transform representation of the 
Green function Foa on Z rf given by: 

, , f exp(2i7rx • u) , v 

^oo(x)= / l \ >-du (96) 

J[0,l] d u l u ) 

where 9(u) = 87 Ylj=i sin 2 (7ru- J ). Remark that is well defined because 
d > 3. We have to prove there exists a constant C > independent of N 
such that 

F N (x) < C(|x| d - 2 + N- 1 ' 2 ) (97) 

Observe that by symmetries of we can restrict our study to the case 
x G [0, iV/2] d . 

We want to show that -Fjv(x) is well approximated by ^(x). We have 
F N ( X ) - i^x) = F N (x) - i£(x) + F^(x) - Foo(x) (98) 

where 

iu , s /" exp(2i7rx • u) 
i£(x) = / -^——^du (99) 

For each k G Z^, we introduce the hypercube Qk = J3^ =1 [k J " / AT, (k J + 
l)/iV) and we divide [0, l] d following the partition U kgZ d Qk- By using this 
partition, we get 

g2«7rk-x/A r g2i7ru-x 



n J «k N (k/N) 

/ 1 1 
f r/np 2l7ru ' x I 



(100) 



Remark that 

2i?rk-x/Af 

due 2 *™'* = — ^^(x/iV) (101) 
^ N 

where 

srn(7ru J ) 



/ 



(102) 



/ \ TT 2i7ru^ TT J11 H""- 
^ U )=ll e 11-^) 

j=l 3=1 V 7 

It follows that the first term on the right hand side of (llOOp is equal to 

(l-<p(x/N))F N (x) (103) 

so that 

Fn (x) = + wmJz L due2m (wo " 

(104) 

The next step consists to show that the second term on the right hand side 
of (jl04H is small. In the sequel, C is a positive constant independent of N 
but which can change from line to line. For each u G Qk, we have 

sin 2 (W) - sin 2 (vrkViV) = vr sin(27r Ci )(u J ' - k j /N) (105) 
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for some Cj G [k J /N, (k J + 1)/N). Consequently, we have 

C 

| sin 2 (vru J ') - sin 2 (irW /N) \ < — | sm(-irk j /N)\ (106) 

Moreover, there exists a positive constant C such that 

Vk G Z%,Vu e Qk, 9 N (u) > C0 N (k/N) (107) 

It follows that the modulus of the second term on the right hand side of 
(I104j) is bounded by 

C * 1 v N-i\sm(nV/N)\ 

l^ml^ — r„(\cIn\2 — (108) 



Since the modulus of the function f(u) is bounded below by a positive 
constant on [0, l/2] d ,this last term is of the same order as 

'sin(7ru : ')| 

'-du (109) 



r- 

Elementary standard analysis shows that this term is of the same order as 

rl yd 

J (N^+r 2 ) 2 

For d > 4, this term is clearly of order iV _1 . For d = 3, the change of 
variables r = N~ l l 2 v gives an integral of order N~ l \ogN . In conclusion, 
we proved 

Moreover, it is not difficult to show that 

|^(x)-F£(x)|<CW-V2 (112) 

Since we have (cf. [14j . theorem 4.5) 

F^x) <C\x\ 2 - d (113) 

we obtained the first point of the lemma. 

For the 1 and 2-dimensional estimates, we have that |Gat(x)| < Gn(0) 
and by standard analysis, there exists a constant C > independent of N 
such that 

Gjv(0) <C f dk (114) 

-Wl" iV-i + E,tism 2 (7rk) 

By using the inequality sin 2 (7rn) > 4-u 2 , one gets G/v(0) is of same order as 

f dk 1 1 - - (115) 



This last quantity is of order y N if d = 1 and log iV if d = 2. 
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Let us now prove the final statement. Assume d > 2 (the case d = 1 can 
be proved in a similar way). We have 



\G N (x + ei) + Gtv(x - ei) - 2Gjv(x)| 

9 e -2«7rk-x/7V 

^E ke zi(l + e 2i7rkl/JV )sin 2 (7rk 1 /iV) 



(k/JV) 



4 

< (27)- 1 



sin 2 (vrki/Af) 



(116) 



□ 



7. Appendix: Equivalence of ensembles 

In this part, we establish a result of equivalence of ensembles for the micro- 
canonical measure < • > N ,s since it does not seem to appear in the literature. 
The decomposition in normal modes permits to obtain easily the results we 
need from the classical equivalence of ensemble for the uniform measure on 
the sphere. This last result proved in [10] says that the expectation of a 
local function in the micro canonical ensemble (the uniform measure on the 
sphere of radius \f~k~ in this context) is equal to the expectation of the same 
function in the canonical ensemble (the standard gaussian measure on M°°) 
with an error of order k^ 1 . In fact, the equivalence of ensembles of Diaconis 
and Freedman is expressed in terms of a very precise estimate of variation 
distance between the microcanonical ensemble and the canonical ensemble. 
In this paper, we need to consider equivalence of ensembles for unbounded 
functions and to be self-contained we prove in the following lemma a slight 
modification of estimates of [10J. 

Lemma 4. Let X n ^ a be the uniform measure on the sphere 

S? nl/2 = j(xi, . . . )Xn ) € R n ; x i = 

of radius r and dimension n — 1 and A£° the Gaussian product measure with 
mean and variance r 2 . Let 9 > and <p a function on M. k such that 

(k \ ® 

X>|J , C>0 (117) 

There exists a constant C (depending on C,8,k,r) such that 

hmsup n \\ n rn i,M - \?{<P)\ < C (118) 

n— >oo 

Proof. This lemma is proved in [10] for <f> positive bounded by 1. Without 
loss of generality, we can assume r = 1 and we simplify the notations by 
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denoting A n nl/2 with A™ and A£° with A°°. The law of {x\ + . . . + Xk) 2 under 
A n is n times a 0[k/2, {n — k)/2] distribution and has density (cf [10J) 



1 r(n/2) (u\ {k/2) - X ( 

f (u) - l {0 <„< n} • - T(k/m{n _ k)/2] y 



(fc/2)-l , ^ ((n-fc)/2)-l 

1 - 

Tlj 



(119) 

On the other hand, the law of (xi + . . . + Xk) 2 under A°° is xt with density 

(cf m) 

g(u) = h/o \ - e -u/2 u {k/2)-l (120) 

yy ' 2 k / 2 r(k/2) y 1 

With these notations, we have 

/>oo 

|A"(0) - A 00 (0)1 < C / u e \f(u) - g(u)\du (121) 
J o 

The RHS of the inequality above is equal to 

f°° / f(u) \ + f°° 
2C u 9 g(u)du + C u e {g{u) - f(u))du (122) 

Jo \SW J Jo 

(f (w) \ 
— — — 1 ] < 2(k + 3)/(n — fc — 3) as soon as fc G 
SW / 

{1, . . . , n — 4}. The second term of (|122p can be computed explicitely and 
is equal to 



r((20 + fc)/2) 

r(fc/2) 



n°T(n/2) 



r (9 + „/2)_ (123 > 

A Taylor expansion shows that this term is bounded by C jn for n large 
enough. □ 

We recall here the following well known properties of the uniform measure 
on the sphere. 

Lemma 5. (Symmetry properties of the uniform measure on the sphere) 
Let Aj? be the uniform measure on the sphere 



ok 



|(x lj ...,x fc )€(M d ) fc ;5;x2 = r 2 | 



of radius r and dimension dk — 1 . 

i) A^ is invariant by any permutation of coordinates. 

ii) Conditionaly to {xi, . . . , x^.}\{xj}, the law of Xj has an even density 
w.r.t. the Lebesgue measure on M. d . 

In the same spirit, we have the following lemma. 

Lemma 6. Let be the uniform measure on the surface defined by 

{k k \ 

(xi,...,x fc )€0R d ) fc ; 5>, 2 = r 2 ; = I 
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We have the following properties: 

i) irfj: is invariant by any permutation of the coordinates. 

ii) If i ^ j G {1, . . . , d} then for every h, I € {1, . . . ,k} (distincts or not), 

Hi) If h / £ £ {1, . . . , k} and i £ {1, . . . , d}, 



fj, r (x l h xi) 



dk(k - 1) 



k - 1 



(124) 



Lemma 7. (Equivalence of ensembles) 

Consider the (a, z/) -harmonic case. There exists a positive constant C 
C(d,£) such that: 



a) ifi + j, 



PoPe! 



d 2 



< 



(PoPejPoPei 



< 



c 

C 



Hi) For any i and any y G Zjy ; we Ziaue 
q^qie.-qiJCp'o) 2 



< 



ClogiV 
N d 



Proof. Let us treat only the unpinned case v = 0. The pinned case is simi- 
lar. We take the Fourier transform of the positions and of the momentums 
(defined by (|70|)) and we define 



q(0 = (i - <s(OMOq(0, p(0 = ^ d/2 (i - <*(0)p(0, 



(125) 



where cj(£) = 2N~ d / 2 \J a Ylt=i sfn 2 (7r£ fc /iV) is the normalized dispersion 
relation. The factor 1 — <5 in the definition above is due to the condition 
^ x p x = ^ x q x = assumed in the microcanonical state. Then the energy 
can be written as 

Kn = ^£{lp(OI 2 + lq(OI 2 } 

= \ E {^ 2 (P(0) + W(P(0) + ^e 2 (q(0) + 3m 2 (q(0} 

Since p x ,q x are real, !SHe(p), !SHe(q) are even and 3m(p), 2Jm(q) are odd: 

9te(p)(0 = «e(p)(-0, J«e(q)(0 = *e(q)(-0 
3m(p)(0 = -Jm(p)(-£), 3m(q)(£) = -3m(q)(-f) 



(126) 



On Z^\{0}, we define the relation £ ~ £' if and only if £ = Let U d 



iV 



be a class of representants for ~ (U^ is of cardinal (N d — l)/2). With these 
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notations and by using (|126|) . we have 

n N =Y. l^ e2 (p(0) + 3m 2 (p(0) + ^ e2 (q(0) + ^ 2 (q(0} (127) 

It follows that in the microcanonical state, the random variables 
((9tep)(0, Pmp)(0, (IHeq)^), Jmq)(0) $6U , 

are distributed according to the uniform measure on the sphere of radius 
V N d £ (which is not true without the restriction on the set Ujy). The classi- 
cal results of equivalence of ensembles for the uniform measure on the sphere 
(|10j) can be applied for these random variables. 

i) By using inverse Fourier transform and (1126j) . we have 

2iirei ■ (77 + 77') 



A 1 



PoPeJ ) =jm E (pW^pW^OL^ n 



N 2d 

N ' £ ?,£',»7,f//0 



(128) 

It is easy to check by using (ii) of lemma [5] that the only terms in this sum 
which are nonzero are only for £' = — £ and r] = —rf. One gets hence 

(PoPe,) 2 ) =^E(|p J (0| 2 |p i W|\ £ (129) 



Classical equivalence of ensembles estimates of |10j show that this last sum 
is equal to {£/df + 0{N- d ). 
ii) Similary, one has 

(PoP^PoPej) 



1 / 9?7rei \ 

= Jp d E (p\0mW(v)&(v')) N , E e W (_L (£<+,/)) 

(130) 

It is easy to check by using (ii) of lemma [5] that the only terms in this sum 
which are nonzero are for £' = — £ and rj' = —77. One gets hence 

iPoPUP^L = J^S E (|P 4 ^| 2 I^Wll/Xp (* + !/)) 

(131) 

Using classical equivalence of ensembles estimates ([TO]), one obtains 

£2 I 1 2 

N I +0(iV^) = 0(iV- d ) (132) 



(PoP^Pi,)^ = ^ 
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iii) By using the symmetry properties, we have 

<^(0?(Op%)iW)>„ >c = o 

for £ ^ — £' or r] ^ —r]'. Hence one has 



^Vo) 2 



1 V" I ~i(c\~uci\~i( \~U i\\ exp (-2i7r(£ • z + £' • y)/iV) 



1 /l ~3(t\~3(n >i 2 \ exp (-2f< • (z - y)/iV) 



Estimates of JTQ] give 

((^(O)V(ei)) 2 ) 



£ 

'> £ V d 



< 



C_ 

~N~d 



It follows that 

q J y (qi ei -0(P J o) 2 
where 



£2 e -2i7re-(-ei-y)/JV _ g-2Mr£-(ei-y)/W 



E 



&0 



+R N 



\Rn\< CN- 2d Y^ 



Edn(27r£7JV)| 



^4aEtisin 2 KVA r ) 



To obtain iii) observe that 
1 

Jq2d 



1 e -2i7rf(-ei-y)/iV _ g-2i7r£-(ei-y)/iV 

Z> ^2 = Tjv(y) 



and 



^4«EtiSm 2 KVA r ) 



log N/N, d = 1 
1/N, d = 2 
l/N d , d>3 



□ 
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